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Abstract 

An algorithm for calculating two-loop propagator type Feynman diagrams 
with arbitrary masses and external momentum is proposed. Recurrence relations 
allowing to express any scalar integral in terms of basic integrals are given. A 
minimal set consisting of 15 essentially two-loop and 15 products of one-loop basic 
integrals is found. Tensor integrals and integrals with irreducible numerators are 
represented as a combination of scalar ones with a higher space-time dimension 
which are reduced to the basic set by using the generalized recurrence relations 
proposed in |Q] . 
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1 Introduction 



Mass effects play an important role in confronting experimental data obtained at the 
high-energy colliders, like LEP and SLC, with theoretical predictions. Precise de- 
terminations of many physical parameters in the Standard Model (SM) require the 
evaluation of the mass dependent radiative corrections. 

At the one-loop level the problem of calculating Feynman integrals with several 
masses was in principle solved in . The history of two- loop calculations is rather long. 
First results were reported more than fourty years ago [[| , |3J . During the last years a 
number of approaches for the evaluation of two-loop diagrams with massive particles 
were proposed || - ||12|| . Integrals with specific topologies or mass combinations were 
considered in |13[ - [p6 |. Most results for scalar master integrals were obtained using 
dispersion relations. The first generalization of the method of integration by parts 
[f27jj to massive integrals was done for the on-shell diagrams with one massive particle 



2q| . The evaluation of the two-loop correction for the photon propagator with the 



aid of recurrence relations was performed in ||17|| . Up to now all attempts to extend 
the standard method of integration by parts |27j to two-loop diagrams with arbitrary 
external momentum and all masses different were unsuccessful. For a short review of 
some techniques for calculating two-loop diagrams, see for instance Ref. . 

Despite these efforts even for propagator type two-loop diagrams no satisfactory 
solution similar to the one-loop case is available so far. In the SM, due to the compli- 
cated structure of integrals and the large number of diagrams (typically thousands), 
no complete calculation of a two- loop self-energy has been carried out. Many differ- 
ent species of particles with different masses have to be taken into account and this 
makes the evaluation of two-loop integrals a rather difficult problem. Existing numer- 
ical methods for evaluating two-loop diagrams cannot guarantee required accuracy for 
the sum of thousands of diagrams. More promising may be the semi-analytic method 
proposed in Ref. ||. It is based on momentum expansion, conformal mapping and the 
use of Pade approximants. 

In calculating Feynman diagrams mainly three difficulties arise: tensor decom- 
position of integrals, reduction of scalar integrals to several basic integrals and the 
evaluation of scalar basic integrals. 

The first two problems were addressed in Refs. ||, |10[ where an algorithm for 



the analytic calculation of two-loop diagrams with several masses was proposed. This 
algorithm, in principle, allows one to express any diagram in terms of several scalar 
integrals but in our opinion the problem was not solved completely. In calculating the 
tensor integrals, the authors of Ref. have to introduce in addition to the original five 
propagators with different masses a fictitious massless propagator. As a consequence, 
in the final result there are not only integrals with the original masses but in general 
also integrals having an extra massless propagator. These integrals may correspond 
to diagrams having a topology which is not inherent to the original diagram under 
consideration. The minimal set of independent basic integrals was not found in M. 
In fact, the algorithm of Ref.|| allows only partial reduction and reveals even in the 
simplest case with one massive particle an overcomplicated result. For example, the 
two- loop correction to the photon propagator in || is given in terms of 4 two-loop scalar 



integrals, though there must be only two [17]. The position of the nearest threshold 



singularity of the two additional integrals is not characteristic to the two-loop photon 
propagator. Without further simplification of the result additional integrals may cause 
potentially large numerical cancellations. 

Recently in Ref. |HJ a framework for treating all two-loop diagrams which occur 
in a renormalizable theory was developed. To reduce tensor integrals to scalar ones 
the authors introduced as a basis ten scalar integrals. For nonrenormalizable theo- 
ries they predict the impossibility to isolate a finite number of functions from which 
all other can be obtained by differentiation. As will be seen from our paper there is 
no difference between two-loop propagator integrals occurring in renormalizable and 
non- renormalizable theories. Any two- loop, so-called "London transport" (or "sun- 
set" ) type diagram with arbitrary masses considered in |30| can be expressed as a 
combination of four two-loop integrals and products of one-loop tadpole integrals. 

In the present paper we propose a method which completely solves for propagator 
diagrams the first of the two aforementioned problems. We formulate an algorithm 
for transforming tensor integrals into a combination of scalar ones. In fact, this is a 
specification to the two-loop case of the general algorithm given in Ref. |]J . We show 
how to reduce all scalar integrals occurring in the computation of an arbitrary two- 
loop diagram to a sum over a minimal set of basic integrals with rational coefficients 
depending on masses, momenta and the dimensionality of space-time d. All recurrence 
relations needed for the reduction of a variety of Feynman integrals to basic integrals 
are presented. The minimal set of basic integrals is given. 

The paper is organized as follows. In Sec. 2 we describe an optimal procedure for the 
simplification of the integrand of a diagram. In Sec. 3 a method for the representation 
of tensor integrals and integrals with irreducible numerators in terms of scalar ones with 
shifted d is presented. In Sec. 4 generalized recurrence relations which are necessary 
for the reduction of two-loop integrals to basic ones are given. In Appendix a complete 
set of basic integrals is given. 



2 General notation and integrand simplification 

The subject of our consideration will be two-loop two-point dimensionally regulated 
diagrams with arbitrary masses. In principle, with the method presented in this article 
one can treat the diagram in tensor form. In Sec. 3 an algorithm for the tensor decom- 
position of individual integrals will be given. Technically the evaluation of diagrams in 
tensor form is more involved than the evaluation of scalar diagrams. In practice one 
need to know the coefficient of the particular tensor structure of the diagram. After 
contracting a diagram with an appropriate projection operator, performing traces we 
commonly encounter integrals of the form: 

/ ( <v) = - d [[ m, k i k ^ k ^ w. 

TT J J C x C 2 C 3 C 4 C 5 

where N(kf, fcf, kiq, k 2 q, k\k 2 ) is a polynomial and 




m\ + ie, 
m\ + ie, 



C3 = (fa - q) 2 -mj + ie, 
Q = (fa - q) 2 -m\ + ie. 



c 5 = (ki - k 2 ) 2 



— m\ + ie 



For brevity, we shall omit the "causal" ie's below. Individual integrals can be specified 
by the dimensionality of the space-time, d and by the powers of denominators z/j, 
called indices of the lines. Figure 1 shows a generic topology of the two-loop two- 
point diagram with lines labeled as indicated. We assume that to the ith line here, 
corresponds a factor l/c^\ In the case of integer i/'s any two- loop diagram with different 
topology can also be reduced to integrals of the form ([]]), if we use the partial fraction 
decomposition for denominators having the same momenta but different masses. 



The evaluation of the integrals (|l|) will be performed in several steps. First, we simplify 
the integrand as much as possible, by bringing it to a standard form suitable for further 
calculations. Second, one has to get rid off the numerator N{k\, k\, faq, faq, fafa) 
by introducing scalar integrals with a higher space-time dimension whenever needed. 
Third, by turns, two-loop integrals with five, four and three lines having different 
powers of propagators will be reduced to a set of basic integrals. Fourth, integrals in 
higher dimensions must be reduced to a basic set of integrals in the generic dimension 
d. Finally, one-loop propagator and tadpole integrals should be evaluated. 

Putting the integrand into a useful form is just a matter of tedious algebra. We 
start the simplification with the repeated use of the substitutions: 




Fig. 1. The master two 



loop two — point diagram 




(k 2 ) 



(k 2 r ( c 2 + m 2 




(*?) 



C 2 \ ^2 

(k 2 ) a ( Cl + m\ 



2 \ <a,u\> 



(2) 



c? [ k\ 



where we used the notation 



< a, Vi >= mm(o;, Vi). 



After performing these substitutions scalar products in the numerator remain only in 
the case when at least one line is eliminated. For integrals of this kind the simplification 
of the integrand proceeds as follows. 

The scalar product (k\k 2 ) remains in the numerator only if invariant C5 is cancelled, 
i.e. the integral is in fact the product of one-loop tensor integrals. In this case the 
substitution: 

tA = i(wj t MM, (3) 

T 

where 

A(h, k 2 ) = k llx (g^ - -^r^j k 2 u, (4) 

allows one to transform the integral into a sum of products of one-loop scalar integrals. 
Integrals with A l (ki,k 2 ), where / is odd (i.e. with odd powers of k\ or k 2 in the 
numerator) are zero since the transverse tensor standing in braces of Eq.(^) will be 
always multiplied by an external momentum q. Integrals with even powers of A(ki, k 2 ) 
are: 

a , ■> r(n + !)rfe 

dVhhik^q^Uk^q^A^^h)- 



r(|)r(n + 4=i) 

x / / d d k 1 d d k 2 f 1 (k 1 ,q,m l )f 2 (k 2 ,q,m l )A n (k 1 ,k 1 )A n (k 2 ,k 2 ). (5) 



Once the factor A(ki,k 2 ) is removed the integrations with respect to k\ and k 2 are 
completely decoupled. 

The integrand can be further simplified by applying the substitutions: 

(ktT _ (*?)« fcs + 2k lQ - g 2 + m|\ <a >" 3> 



c 3 c 3 \ K l / 

(kl) a _ {klT { c, + 2k 2 q-q 2 + ml \ <a ' Vi> 
C4 4 C4 4 \ k 2 J 



(7) 



After these substitutions have been made, kf, k\ will remain in the numerator only 
in bubble-like integrals. Integrals with irreducible numerators reveal themselves as 
integrals having scalar products {k\q) a and/or {k 2 qY in the numerator. In principle 
it is possible to eliminate {kiq), (k 2 q) at the expense of kf, k\. As we shall see later, 
integrals with (kiq), (k 2 q) will produce scalar integrals with a smaller shift of d (and 
therefore simpler expressions for further computation) than integrals with kf, k\. 

If C3 and C4 are absent in the integrand then the scalar products (kiq), (k 2 q) can be 
eliminated. In this case, for an arbitrary scalar function f(ki, k 2 ), 

d d k 1 d d k 2 f{k 1 ,k 2 ){k iq yi(k 2q y> = v(u u v 2 ) t (8) 

is zero if v\ + v 2 is odd. For V\ + v 2 even, the scalar products (kiq), (k 2 q) can be 
eliminated from the integrand applying the recurrence relations derived in [[3TJ : 



(d+ v x + v 2 - 2)v(v x , v 2 ) = g 2 {(z/x - 1)^1'+ v 2 (k x k 2 )2-}\- o v(v x , v 2 ), (v x > 1) (9) 



{d+is 1 + v 2 -2)v(v 1 ,is 2 ) = q 2 {(u 2 -l)k 2 2 2- + u 1 (k 1 k 2 )l-}2- ov^,^), {y 2 > 1), (10) 

where l^v^i, u 2 ) = f (^1 ± 1, ^2) etc., and the o sign means that factors kf, k 2 , k\k 2 in 
braces must be considered under the integral sign in v. To simplify one-loop tadpole- 
like integrals, rather frequently, application of the simple formula 

turns out to be useful. 

Using the formulas presented in this section, a scalar contribution ([!]) will be trans- 
formed into a sum of integrals with irreducible numerators {kiq) a {k 2 q) 13 having at 
maximum four lines and a variety of scalar integrals with different indices. Integrals 
with irreducible numerators (kiq), {k 2 q) will be considered in the next section. Notice 
that, the integrals (§) can also be regarded as integrals with irreducible numerators 
and they can be treated according to the formalism discussed in the next section. 



3 Tensor integrals and irreducible numerators 

Tensor integrals can be written as a combination of scalar integrals with shifted space- 
time dimension multiplied by tensor structures made from external momenta and the 



metric tensor. This was shown in Ref. [32 for the one-loop case and in Ref. M for an 



arbitrary case. In the approach of Ref. |jj the reduction of tensor integrals to such a 
representation is performed by applying a tensor operator to a scalar integral. For the 
two-loop tensor integrals of interest the following relation holds: 

d d hd d k 2 , , , , m . rol rr d d k l( i d k 2 



"\ , "< <",'% ^iM • • • • • • &2A S — ^jti...A s (<?, {dj}, d + ) / / , ,',»'.-. ' 

c l c 2 c 3 c 4 c 5 J J c l c 2 c 3 c 4 c 5 



(12) 



where 



A 9 



dm 2 ' 



and d + is the operator shifting the value of the space-time dimension of the integral 
by two-units: d + /^ = l( d+2 \ On the right-hand side of Eq.(|l2|) it is assumed that, 
at the beginning, invariants Cj have different nonzero masses and after differentiation 
with respect to mf they must be set to their original values. 

To illustrate the method of Ref. [|IJ we will derive here an explicit expression for 
the tensor operator I) tl ...A s (?, {dj}, d + ). The main ingredients of the derivation are 
independent auxiliary vectors 01,02 an d the use of the a- parametric representation 
[Q . The tensor structure of the integrand on the left-hand side of fll2|) can be written 
as 



Id d d d 

fci w • • • kin r k 2Xl . . . k 2Xa = — -= . . . -7; ■= . . . ^ exp [i(oifei + 02^2)] 

t r+s da lfll da^ r da 2Xl da 2Xs 



ai=Q 

(13) 



To convert the integral 



G {d \q 2 ) = I I Jl^tl V5 exp [i( ai kt + a 2 k 2 )] . (14) 



Ci c 2 c 3 c 4 c 5 



33|1 ). Transforming all propagators into a parametric form 

1 9— v roo 



into the o-parametric representation we apply standard methods (see for example, 

ametric form 

da a" -1 exp ia(k 2 — m 2 + ie) , (15) 

(16) 



(k 2 -m 2 + ie) u T(u) Jo 
and using the d- dimensional Gaussian integration formula 



J d d kexp \i(Ak 2 + 2(pk)) 



d 


• 2 " 

IP 


=i (w exp 


~~A 



we can easily evaluate the integrals over loop momenta. The final result is: 

Q(a, at, o 2 ) 



jv\ d JL i "> r°° r°° docjOi? 1 



/ / r(^)Jo Jo [D(aj] 



-exp 



m l —it) 



(17) 



where 



-D(o) = o 5 (oi + o 2 + o 3 + o 4 ) + (oti + o 3 )(o 2 + a 4 ), 



Q(a, ai, a 2 ) = [(ai + o 2 )(o 3 + 0:4)0:5 + oio 2 (o 3 + o 4 ) + o 3 o 4 (oi + o 2 )]g^ 



+ (qai)Qi + (ga 2 )Q 2 + a\Q n + a 2 Q 22 + {a x a 2 )Q 



12; 



(19) 



and 



Qi 

-4Qn 
-4Q22 
-2Q12 



03O5 + 04O5 + o 2 o 3 + 0304, 
04O5 + 03O5 + 01O4 + 03O4, 
o 2 + o 4 + o 5 , 
Oi + o 3 + o 5 , 

05. 



(20) 



It is straightforward to work out from (|I2"D , ( |i3|) and the explicit expression (]ITD that 

9 



1 r d 

n 



n 



2 An 



x exp 



ga^Qi + (qa 2 )Q2 + a\Q n + AQ22 + (aia 2 )(5i 2 ) P 



a J= 



.(21) 



This operator is a particular case of the more general one derived in Ref. Q. Formula 
flTjp can be used for the direct evaluation of two-loop tensor integrals. Notice that in 
order to obtain the tensor decomposition of the integral no contractions with external 
momenta and the metric tensor and no solution of a linear system of equations are 
needed. 

Integrals with irreducible numerators 

d d ktd d k 2 



«' - // 



c l c 2 c 3 c 4 c 5 



^(kiq) r (k 2 qY 



22 



can be regarded as a contraction of the tensor integral (|T^) with . . . q Vs . For the 
scalar integrals (p2]) a somewhat simpler formula, analogous to (|12"D , can be derived by 
introducing auxiliary scalar parameters (3j. Similar to ([13]) we write 



(hq) r (k 2q y 



QT 



s 



exp{i[/?i(/cig) + (3 2 (k 2 q)]} |ft=o 



(23) 



{idfrY (id/3 2 ) 

Carrying through the steps which were leading us to (pT|) , we find the relation: 

ijf =T„( g> {5 i } I d + )J« ) , (24) 

where 

1 <9 r 9 s 



i r+s d(3{ dpi 



exp {ig 2 [Qi/3i + Q 2 /3 2 + Qn/3 2 + Q 22 (3 2 2 + 



(25) 



with Qj, Qr, given in (|20|). 

The evaluation of the scalar integrals in the form fl24"[) is more efficient than the 
use of its tensor analog (|T2"D. Notice that, as we mentioned in the previous section, 
irreducible numerators can appear only if at least one of lines 1 to 4 is contracted. In 
this case the operator T rs (q, {dj}, d + ) becomes simpler because one can set equal to 
zero those a parameters which correspond to contracted lines. One can easily build 
up transformation operators ^...^(g, {dj} or T rs (q, {dj},d + ) by using any computer 
algebra system. Running FORM |34] on a PC Pentium 90 it takes usually several 
seconds (or minutes in complicated cases) to construct these operators. 

As an application of the formalism presented in this section, let us consider a typical 
integral with an irreducible numerator: 



in 



d d k 1 d d k 2 

C 2 CzC^ 



(hq)(k 2 q). 



(26) 



For this particular case, dropping irrelevant terms with f3f,f3 2 in (p5|), the transforma- 
tion operator reads: 



11 



d d 



exp [iq 2 [ftQj + (3 2 Q 2 + /3i/3 2 Qi 2 ] p} 



0i=O 

TT 



(27) 



Since C\ and C4 are absent in the integrand (i.e. corresponding lines are contracted), 
we have to set a± = 0:4 = and thus obtain from (f2"U|): 



Qi = a> 3 (a 2 + a 5 ), 
Substituting these expressions into 



a 3 a 5 , 



Q12 = ~2 a5 ' 



11 



2tt 



we get 
q 4 



d + d 5 + ^-(d + ) 2 d 2 3 d 5 (d 2 + d 5 ). 



7T J 



(28) 



(29) 



With this operator fl2l|) leads to the desired relation: 

(faq)(k 2 q) 



J jd d fad d k 2 - 



C2C3C5 

rd d+2 fad d+2 k 



2tt 2 



c 2 c 3 cl 



2 + q r 

TV 



J Jd d+ %d d+ % 



C2 ^3 ^2 
2 C 3 C 5 



+ 



c 2 44 



(30) 



Integrals on the right-hand side of (|30|) can be reduced to basic ones in the generic 
dimension d by using recurrence relations given in the next section. Notice that for 
these integrals at arbitrary d an analytic expression in terms of Lauricella functions is 
available |19[ and therefore formula (|23|) shows that the "sunset" type diagrams with 



any number of scalar products in the numerator are always expressible in terms of 
restricted number of the Lauricella functions. 



4 Recurrence relations 

In the present section we will be concerned with generalized recurrence relations. Ap- 
plying methods presented in previous sections we have reduced the problem of calcu- 
lating two-loop diagrams to that of evaluating scalar integrals without numerator and 
having different shifts of d. The reduction of scalar integrals with different powers of 
propagators (or different indices) to a sum over a minimal set of basic integrals in an 
arbitrary dimension will be done by means of generalized recurrence relations. The 
method of their derivation was described in Ref. For the two-loop general mass 
case the derivation of these relations is rather tedious and for brevity of the presen- 
tation will be omitted in the present paper. We will give only the final formulae and 
describe their application. To avoid overcomplication by irrelevant indices it is conve- 
nient to introduce a separate notation for two-loop integrals with five, four and three 
propagators: 

F (d) J_ f f d d fad d k 2 



[(fa - q) 2 - ml] U3 [(k 2 - q) 2 - m%\"*[(ki - k 2 ) 2 - m 2 ]" 5 
d d fad d k 2 



"> J J (/,•,•-/,•,)-' ///-/' 

1 



[hi - ml] y *[(fa - q) 2 - ml]^[(k 2 - q) 2 - ml]"*' 

/(<*) (a 2) = 1_ ff d d fad d k 2 

•V^ 3 W ) nd JJ [£2 _ m 2], x -k 2 ) 2 - ml]"* [(k 2 -q) 2 - m§]«* ' 1 ' 

One-loop integrals will be denoted as 

= — [ (vr\ 

UlU2 niJ [kj-m 2 ]^ [(fa - q) 2 - m 2 2 p ' 1 ' 

T - ) = ^f \ [k\-m\\^ (33) 



Where no confusion can arise, we simply refer to the above functions as F^, V^ d \ . . . 
etc. 

The application of recurrence relations to will produce V^ d \ and more 
simple one-loop integrals. In turn integrals will produce plus one-loop inte- 
grals. Thus one should first apply the recurrence relations to F^ d \ then to and 
after that to j( d \ 

4.1 Integrals with five propagators 

The integrals originate from diagrams with the topology given in Fig. 1. The first 
step will be the reduction of indices v\, . . . , 1/4 by using the recurrence relations which 
derive from the generalization of the method of integration by parts [27[] to massive 
integrals. It is particularly convenient to represent the recurrence relations in terms of 
shift operators decreasing or increasing indices of the integral by one unit |27[|: 



1 ± p{d) F (d) f (OA) 

with a similar convention for V^ d \ and other integrals. Notice that an analogous 
notation we already employed in Section 2. 

The indices ui, . . . , 1/4 can be reduced to or 1 by iteratively applying the relation: 

2v l Al + FW ViVB = Ud - 2v x -u 3 - z/ 5 )A 1 + A 345 h5 + (2- - 1~) - z/ 3 3+:r] 



+ A 2 [^l + (5- - 2") + z/ 3 3 + (5~ - 4~) + u 5 - v x ] 

+A 6 hl + 3- + z/ 5 5 + (3- - 4-) + v 3 - ^]} F^ V2l/3W (35) 
along with relations for 2 ± F^\ 3 ± F( d \ 4 ± F^^ following from ( |35l) by replacements: 



2+F 
4+F 



(ui, V3, mi,m 3 , 1 ± , 3 ± ) <-> (u 2 , v±, m 2 , m 4 , 2 ± , 4 ±N 
(z/i, z/ 2 , mi, ma, 1 ± , 2 ± ) «-> (z/ 3 , z/ 4 , m 3 , m 4 , 3 ± , 4 
(z/i, V2, mi, m 2 , 1 ± , 2 ± ) <-> (z/ 4 , 1/3, m 4 , m 3 , 4 ± , 3 ± ), (36) 



which are due to the symmetry of the integral. In the above formula we used the 
following abbreviations: 

A 42 17 2 1 2 1 2 1 2 2\2 / 2 2\/ 2 2m 2 

= g m 5 — [(m 1 + m 2 + m 3 + m 4 — m 5 )m 5 — (m 4 — m 3 J(m 2 — m-JJg 

, / 2 2\( 2 2\ 2 i / 2 2 2 , 2\/ 2 2 2 2\ /o^N 

+ (m 3 — m 1 J(m 4 — m 2 )m 5 + [m l — m 2 — m 3 + m 4 ){m 1 m 4 — m 2 m 3 J, (37) 



<9A 



dmj 



2 ' 



A, 



<9A 

dq 2 



o 2 2 2/ 2 , 2 , 2 

2q m 5 — m 5 (m 1 + m 2 + m 3 



+ 777-4 - m 5; 



+ (mi 



m 



;)(m 



777 



A 4 i 4 i 4 o/ 2 2 , 22, 2 2\ 
ijk = m\ + + m k — 2\m i m i + m i m k + m^m k ) 

= (rrii + frij + mj,)(mj + m 3 - — m k )(mi — rrij + mk){rrii — m,- — m&) 

^ijki^jik ^kij) ^jik^kiji 



(3f 



Uijk 



-<);A;, 



ijk 



111 

m\ - - m fe , 



A 



U'6 



4 i 4,4 o/ 22, 22, 2 2^ 

mj+m-+g — /(ifij/n , + m { q + m j q) 



(39) 

It is worthwhile to note that the repeatedly appearing expression (|3"T| ) is the well known 
Cayley kinematical determinant: 



A 



1 

1 

1 q 2 
1 m\ 
1 m\ 



1 1 
m\ 777,3 



m\ 



7777 



mt, 



m\ 



m\ mg 







(40) 



To reduce the index z/5 another recurrence relation is needed. It has a form similar to 
551) 



2^A5 + Fj?) 2! , 3 ^ 5 



{A 136 [z/ 3 3 H 



z/il + (2- -5-)] 



+ (d - v x - is 3 - 2u 5 )A 5 + A 2 [^ 5 5 + (l" - 2~) + z/ 3 3 + l" + v x 



+A 4 M + 3- + ^ 5 5 + (3- - 4-) +1/3- ^]}4t,3^s- ( 41 ) 

If for some mass values A = 0, then the formulas (|35|), (^TJ) are of no use. This case 
needs a separate treatment. For arbitrary q 2 , the relation A = can be fulfilled only 

at 

777^ = 0, m 2 = m^, m\ = m\. (42) 

If this condition is satisfied, substituting ( f4~2"D into (|35"D (or (f41~D ) yields the simpler 
relations: 



(d - v x - u 3 - 2u b )Fl% 



hi + (5" 



^ 3 3 + (5--4-)]F(f) 



^1^2 VZ U4, Vr 1 



(43) 



(d - 1* - 1/4 - 2u 5 )F^ U2U3UiU5 = h2+(5- - 1-) + z/ 4 4+(5- - ^)]F^ mvt . (44) 

By repeated use of these relations, the integrals can be reduced to a combination 
of integrals and one- loop ones. Prior to applying fl43|) , ([44"D , in some cases, it may 
be more efficient to first use the following relation which reduces the sum of indices of 
the integral 

A 136 l + F,ll 



1/11/2^3^4^5 



{u 613 [za1 + 3- + z/ 5 5 + (3--4-)] 



2tj 136 [z/ 2 2 + (5- 



i/ 4 4 + (5~ -3" 



V\ + 3i/ 3 ] 



+2m2[z/ 5 5+(l- - 2-) + z/ 3 3+l- 



^i-^ 3 llFM 



J !/l»/2^3^4^5 ' 



(45) 



Analogous recurrence relations for follow from (g|) and substi- 

tutions (^). Here, and in the sequel, we use the notation 

i 

Ei = 3d - 2 53 Vj ~ 2- (46) 

Also for particular values of the external momentum squared can be reduced to 
simpler integrals. Again this happens if the determinant ( (4"0"| ) vanishes. If m 5 ^ 0, then 
A = is a quadratic equation for q 2 with solutions 



q± = m\ + ml - — | [w 2 i 5 W435 ± \j A 2 i 5 A 435 ] . (47) 
Furthermore, if m\ = 0, m\ ^ m 2 , m% ^ then we have A = at 



o (m 2 — m 2 , — m^ + m 2 )(mlm 2 — m^m 2 ) 

/ 2 2N / 2 2\ ' ( 

[m2 ~ mf){m^ — mj 



At A = the left-hand sides of Eqs. (^), (f4l|) are zero and we are left with recurrence 
relations of a form similar to (f£|), (@)- By repeated application of this kind of relations 
again is reducible to J^ d > and one-loop integrals. For example, the special point 
mi = m 4 = 0, m 2 = m 3 = m 5 = m is the on-shell value of the momentum 



2 2 

q = m . 



This kind of integrals are encountered in calculation of the on-shell diagrams for the 
fermion propagator of QED [p8|j . Simplification of integrals at some specific values 



of q 2 can be useful, for example, in establishing gauge invariance of physical ampli- 
tudes. Considering special cases we have been tacitly assuming that the dimensional 
regularization does not break down at the considered specific values of q 2 . 

Now we proceed with a few remarks concerning a reduction of the space-time di- 
mension d. If we simplify the numerator of a scalar integral according to the rules 
described in Sec. 2 then the integrals F^ d > will occur in the final result without change 
of dimension. Decomposing tensor integrals we will encounter in general F'^'s with 
shifted d. The connection between Fif xll and -Fifm which suffices to reduce these 
integrals to generic dimension can be worked out from the relation derived in Ref . |]J : 

F (d-2) = D(Q\p{d) U§\ 

Here D(d) is a differential operator obtained from ( |1~8| ) by converting <x,- into dj. By 



applying the recurrence relations already presented in this section Eq. (|49| ) can be 
written in the following form: 

A^nrn = -q 2 (d-3)(d-A)Fif ) in + (P 1 l- + P 2 2- + P 3 3- +P 4 4- + P 5 5-)Fifl u , (50) 

where P{ are polynomials in masses, q 2 , d and operators l + ,..,5 + . The polynomial Pi 
does not depend on the operator i + . The final formula connecting F^ n and -Fnm 
with all masses arbitrary is rather long and for this reason it will be not given here. 
The reader can easily obtain it by using the recurrence relations given in the present 



article. Notice that at A = formula (|50|) gives us a relation for -Finn m terms of 
simpler integrals. 

To conclude this subsection, we note that in the case of arbitrary masses and ex- 
ternal momenta by applying (|35|), ([41]) , integrals F„fl 2l , 3U41/5 with integer Uj > will 
be expressed in terms of -Finn, K-il^^' ^u% 2 v3 anc ^ more simple ones. Notice that 
only one integral with five propagators, namely Finn is irreducible with respect to 
more simple ones. It will be the first representative of our minimal basis of two-loop 
integrals. As mentioned previously in the special case of all integrals F^ and V {d) 
are reducible to and one-loop ones by means of (f43|). 

4.2 Integrals with four propagators 

Now we shall consider the integrals V^ d ' which at first sight look simpler than the 
integrals F^ d \ but in fact their reduction causes more problems. The reason is twofold. 
First, the V {d) 's are less symmetrical than the F^ 's. Therefore the number of different 
relations needed is larger. Second, there are more cases to be distinguished because 
different mass combinations lead to different recurrence relations. The topology of the 
diagram corresponding to the integrals V^ d ' and the labeling of lines are shown in Fig. 
2. 




Fig. 2. Two — loop diagram with four propagators 

In many cases the most efficient way to reduce integrals is to first apply the recurrence 
relations which reduce the sum of indices of the integral. Relations of this kind, which 
reduce v 12 to or 1, are 

Ai34^1 + l$U, 4 = {2m^ 3 3 + (l- - 4-) + M4 i3^il + (3- - 4-) 

+u lu (d -vi- 2^ 3 ) + 2m 2 3 (u 1 - u 3 )} V^ m , (51) 



A 246 u 2 2 + vS 2 u,u 4 = {2m^il+(4- - 2") + 2m^3+(4- - 2") 
- 2m^ 4 2^4 + + (2m\ - u i2fi )v 2 2 + 4~ - (E 4 + 2)(4~ - 2~) 

+u 246 (d - 3v 2 ) + 2m^K - ^2)} V^ Vi . (52) 

The index 1/3 is lowered to zero or one by using the relation which follows from ([51]) by 
interchanging 

v\ <-> 1/3, mi <-> 777.3, l 1 * 1 ^ (53) 



An alternative way of reducing v 3 is to use first the relation 

2ml^3 + yM,^ = K 13 z/!l + + ^l + (4- - 3") + d - ^ - 2u 3 ]V u % sU4 , (54) 

and then to apply fl5*T|). 

The reduction of v± is more problematic. It can be done at the expense of increasing 
V\ and v 2 by applying the relation 

2mlv 4 4 + V$ 2V3V4 = {m 3 i4^i1 + + u G24 u 2 2 + 

+z/xl + (3- - 4") - v 2 2+4- + d-u 1 -u 2 - 2z/ 4 } V^ m - (55) 

In turn, u 12 are to be lowered by means of (|5"ID, (|52|). As an alternative we may use 
the more symmetrical relation: 

2mlv 4 4+V^ m = {-2m 2 v 1 l+ + u 6M v 2 2+ 

-2u 3 m 2 3 Z + ~ u 2 2 + 4~ + 2d- 2v x -v 2 - 2v 3 - 2u 4 ) V^^, (56) 

which increases not only z/ 12 but also v 3 . 

For arbitrary q 2 there are two special cases for which some of the above relations 
do not work. Relation (FTD) which reduces v\ and the analogous one reducing v 3 are 
nonapplicable if the masses m 1 ,m 3 ,m 4 are subject to the condition 

A 134 = (rrii + m 3 + Tn 4 )(mi + m 3 — m 4 )(m 1 — m 3 + 777,4) ( m i ~ m 3 ~ m A) — 0. (57) 

Similarly the relation reducing the index z/ 4 cannot be used if 

m 4 = 0. (58) 

These cases require additional investigations. 

When all masses are different from zero but satisfy Ai 34 = we found several 
recurrence relations which allow us to reduce all V^- integrals to simpler ones with at 
least one line contracted. For example, one possible relation of this kind 

[2m\u lu {d -v x - 2v 3 ) + (z/ 3 - l)ul l3 + 4m\m\{v x - is 3 )}V^v 4 = 

[w4i3^3 + 4~(4~ - 1~) - Am\m 2 3 v 3 Z + l' + (Am\m 2 3 + u\ l3 )v 3 Z + 4' 

+ ^413(^3 - 1)1" - (d - 2vi - v 3 + 1)m 4 i 3 3" 

+ ( d _ 2v x - 2v 3 + l)u A13 4-]V$ 2l/3U4 . (59) 
follows from Eq. (O) with A 134 = and using relation 

2mlu 1 l + V u % 3U4 = [u 413 v 3 3+ + z/ 3 3 + (4- -l-) + d-u 3 - 2v l ]V^ vm , (60) 



which is obtained from performing substitutions (|53|). 

Before applying relation fl59|) the indices ^1,2,4 should be lowered to 1 with the help 
of the relations fl52|), (|55l) and (^). Then applying the relation (|5^) as often as needed, 
one line in will be eliminated. Notice that, the term 3 + 4 ~4~ in (|59|) produces 
factors in the numerator which are supposed to be treated by the method proposed in 
Section 3. 

If one of m 2 , m 2 ,m 2 is zero then A 134 = requires two other masses to be equal. 
There are two distinct cases: 

Case 1: at m\ = (m 2 = 0) and m\ = m\ (m\ = m 2 ) any integral with four 
lines can be reduced to a combination of simpler integrals with three lines by using the 
relation: 

(d - 2v x - Vz)V^ m = v 3 3+(l- - 4T)V^ m . (61) 

Case 2: if m\ = and m\ = m\ then A 134 = and the index v 2 should be lowered 
to one by means of (|52|) . Then the relation 

(q 2 - m\){v x - 2v 2 + 2z/ 4 + 2)4 + V$ aV3V4 = {v x (m\ - u B12 )l + + 4m 2 z/ 2 2+ 
+ 2m\v- i Z + + [£ 4 - 2m\v l l + - 2m 2 z/ 3 3 + ]2-4+ 

+^(g 2 - _ £ 4 } V™^, (62) 

is applied until one of the indices ^2,3,4 will be zero. After each iteration it is reasonable 
to reduce v 2 again, in order to avoid a zero on the left-hand side of ( |62"D . 

Let us consider now the case m\ = and Ai 34 7^ 0. Since Ai 34 7^ the relation 
( |51"|) and its analog for 3 + V^ may be used to reduce z/i >3 to 1. Then, instead of (|55|), 
the relation 



(q 2 - m 2 2 )[(mj - m\){d - 2v 2 - 2u 3 + 2z/ 4 ) + 2m 2 (z/ 3 - v x )\V^ V2 

{(m 2 - mD^mlu^ + 2m 2 z/ 3 3 + - E 4 - 2)(2 _ - 4") - 4m 2 2 v 2 2 + 4r] 

+2{q 2 - m 2 2 )[m^ 3 3 + (l- - 4") + m 2 z/ 1 l + (4- - 3-)]} V^ m , (63) 

must be applied until one of the Vi will become zero. After each application of fl63p, it 
is reasonable to reduce to 1, if they have been increased after the use of flB5|) . 

Again the integrals at some special values of q 2 admit a reduction to simpler 
integrals. These special values are determined by the equation 



A 246 = 0, (64) 



which has the simple solution 



q± = (m 2 ± m 4 ) 2 . (65) 

At A 24 6 = the left-hand side of (|52| ) is zero. With the help of equation ( |55|) we 
exclude the term 2~4 + from the right-hand side of ([52]) and obtain a relation which 
can be used to eliminate one line in V^ d \ After each application of this relation one 
should reduce indices ^1,3,4 to 1 if they have been increased by the previous iteration. 



So far, we have considered only relations connecting integrals with the same space- 
time dimension. The whole concept of the method presented in this article would not 
work if we would not able to express integrals with shifted dimension in terms of those in 
a generic dimension. In fact there are several possibilities to reduce the dimensionality 
of an integral. As we shall see, one possibility always exists and it amounts to first 
express the integrals V^*^l 4 in terms of Vini and simpler ones and then to express 

all Vii^ 2k ^ in terms of l^mi pl us simpler ones. 

The other possibility is to reduce simultaneously the indices and the space-time 
dimension. We found several recurrence relations of this kind: 

4gV^ 3 l + 2+3+^t 4 = 

— {Qd — 6u% — 5z/ 2 — 6z/ 3 — 4z/ 4 — Qm\vil + + (w 624 — 4m 2 )z/ 2 2 + 

-6m 3 z/ 3 3 + + 2(w 624 - m> 4 4 + - ^2+4" - 2z/ 4 4 + 2-} K^L^- (66) 



v^(^+i)^i + i + 2 + v 1 sta M4 = 

- {dd — 6^1 — 5z/ 2 — 6z/ 3 — 4z/ 4 — 2(w 246 + 3ml)vil + 
+ («624 - 4m 2 )z/ 2 2 + - 6m 3 i/ 3 3 + + 2(w 624 - m 4 )z/ 4 4 + 

+2za1+(4- - 2-) - z, 2 2 + 4- - 2^4+2'} V^ m , (67) 

4g 2 m^ 1 z, 2 (z, 2 + l)l + 2 + 2 + ^^ 4 = 

- {[(ml - m 2 3 ) (2d - 2v x - u 2 - 2v 3 - 2z/ 4 ) + m|(2i/i + u 2 - 2v 3 )] 

- Mi 34 [2m^il + + u 246 ^ 2 2 + ] + w 3 i 4 [2m 3 z/ 3 3 + + z/ 2 2 + 4~] 

-27^2+ + z/ 4 4 + )(3- - 1~)} V^ m . (68) 

Exploiting the symmetry of V^ d \ two additional relations follow from (]67), (68) and 
the interchangements (|53|). Whenever possible, the above mixed recurrence relations 
should be used first. Their application is more efficient than the step by step reduction 
of all V ^ 's to master integrals in different dimensions and then the reduction of master 
integrals to the generic dimension. 

The key relationship between and V^ux > which is necessary for the d— re- 
currences, follows from the equation (see Ref. [0): 

n ( m 2) = (did* + did 3 + OA + d 2 d 3 + <9 3 <9 4 ) V}$ v (69) 

To simplify the right-hand side of Eq. fl5"9"| ) we use the recurrence relations fl5T p ,(]5^ ) , 
(p5|). After some additional algebra a relatively cumbersome formula is obtained: 

(3d _ io) (d - 3) 2 m^mi = ^ {(3d - 10)A 134 [A 246 + u 246 (2- - 4")] 
-8(d-A) 2 q 2 mlml4-}v}fn 2) 



+ {(d - 3)ml[Amj(d - 4)1+1+ - 8(d - 3)m 2 2 2 + 2 + + (3d - 10) (d - 3)2+] 

+ u 62 ±[-(3d - 10)m 134 (1 + - 3+)+ (d - 2)(d - 3)m;jl+]2+ 

-m^[(d - 4) 2 Ml34 2+ + (3d - 10) M4 i 3 3+] 1+} ATV^ 

- (d - 3)m^ {4(d - 4)m?l + + (d - 2)[u 624 2+ - 2m^4+]} 1+3-^An 

+ ^(3d - 10)n 13 4K 24 2+- 2m^4+](l-3+- 1+3-)^. (70) 

Notice that all terms with on the right-hand side of (ffOf) are accompanied by 

shift operators contracting one of the lines i.e. producing simpler integrals. We kept 
ViiIy for compactness of the formula, instead of specifying particular integrals obtained 
from it by contracting lines. It is worthwhile to note that in the relation ( |70"D when 
A134 = or A 2 46 = 0, the term Vim with uncontracted lines drops out and we obtain 
expression for V^ m in terms of simpler integrals. 

Recall that if ml = 0, all for arbitrary d are reducible to integrals plus 
simpler ones and therefore in this case the problem of reducing d is transferred to 
simpler integrals. 

Concluding this subsection it is worthwhile to mention that in the general mass 
case, after applying appropriate recurrence relations, the only integral which will 
remain is V^i- This will be the second integral in our minimal basis of integrals. 
In fact in calculating diagrams with the topology given in Fig.l. four integrals Vifli 
corresponding to different mass distributions will occur (see Appendix ). All these 
integrals are to be included in the basis. 

4.3 Integrals with three propagators 

The remaining nontrivial integrals (see Fig. 3), which are of interest in this subsec- 
tion, are rather symmetrical though the corresponding recurrence relations are more 
cumbersome than in the previous cases. 




Fig. 3. Diagrams with three lines 



To reduce integrals with three lines to a basic set of integrals two kind of recurrence 
relations are needed. The first one reduces the sum of indices of the integral if at least 



two lines have indices exceeding 1. The relation applicable when v\ > 1 and z/ 2 > 1 
has the form: 



where 



2z/ 1 z/ 2 D 123 l + 2 + 4fl 2U3 (q 2 ) = [2u 1 h 123 l + + 2u 2 h 213 2 + + 4z/ 3 m 2 cx 123 3+ 
+ v 2 v 3 m 2 3 (j) 2V3 l~2 + 2 l + + z/iZ/ 3 m 2 0i 23 l+2~3+ - 2z/iZ/ 2 pi 23 l+2+3~ 

4%v 3 (q 2 ), (71) 



-~£ 3 (d - I/i - I/ 2 - f 3) 0321 



Ajfc = g 8 - 4g 6 (m 2 + mj + m 2 ) + g 4 [6(m 4 + m 4 + m 4 ) 

+ A(m 2 m 2 + m 2 k m 2 + m 2 k m 2 )} — 4g 2 [m| + + — m 2 (m 4 + m 4 .) 

- m 2 (m 4 + m 4 ) - m 2 (m 4 + m*) + 10m 2 m 2 m 2 ] + A 2 ,,, (72) 



^^ijk 6 1 o 4/ 2 1 2 1 2 

pijfc = -7-^- = -9 + 3g (m, + m i + m fc 



4 dg 2 

g 2 [3(m 4 + m 4 + m 4 .) + 2(m 2 m 2 + m 2 m^ + m 2 m 2 )] + + m® + m| 
m 2 (m 4 + m 4 ,)— m 2 (m 4 + ml)—m 2 k (mf + m*) + l0m 2 m 2 ml, (73) 



_l d ( d | g | g i *M r 

i,? 2 9m 2 \ (9m 2 9m 2 <9m| <9g 2 J 13 
4[g 4 + 2g 2 (m- - m 2 - m 2 ,) + (m 2 - m 2 ) 2 + m 2 (2m 2 + 2m 2 - 3m 2 )] 
= 4[g 4 + 2q 2 u ijk + u ijk u jik + u ijk u kij - u jik u kij ], (74) 

= --(d - i/j - 2vj)<f) ijk - -(d - 2v t - Vj)<f) jik 

--(2d — 2vi-2vj-v k -l)<f>kij, (75) 

h ijk = -~(d - 2vj - v k )m 2 k (j) ijk - -(2d - ^ - 2z/,- - 2v k - l)m 2 (fi kij 

+ (d- Uj - 2v k )p ijk . (76) 
Two more recurrence relations follow from (|7l|) by the interchanges 

l+3+J (d) : I / 2 ,m 2 ,2 ± ^z/ 3 ,m 3 ,3 ± , 

2+3+J (d) : ^i,m 1 ,l ± «-> z/ 3 ,m 3 ,3+. (77) 

The above recurrence relations are of no use if two lines have an index equal to one. 
In this case the second set of recurrence relation is needed. To reduce the index of the 
first line the following relation applies: 



2m 2 D 123 z/ 1 (z/ 1 + 1)1+1+ Jl% u ,(q 2 ) = {"S 3 (d -v x -v 2 - v 3 )p 



123 



+m\m\(j) l23 v 2 v 3 l 2 + 3 + + m^m302i3 z/ i z/ 3l +2 Z + + m\m\(j) 3X2 v x v 2 l + 2 + Z 

+ (d - 2- 2z/ 1 ) J D 123 z/ 1 l + + m^ 123 z/ 1 l + + m^ 213 ^ 2 2 + + m 2 S 312 z/ 3 3 + } 4% 2U3 (q 2 ), (78) 
where 

= -{d - 2x/j - v k )m\<\)j ik -(d- Vj - 2v k )m 2 (f) kij 

+ 2{2d-v i -2v i -2v k -\) Pijk . (79) 

Recurrence relations for the reduction of z/ 2 3 can be easily deduced from (|78|) by sym- 
metrical interchanges similar to fl77|). 

We see from (|78|) that if two indices are equal to 1 then it is impossible to reduce 
the third index to 1. This means that our minimal set of integrals will include not 
only J X1X but also J 2XX , j[ 2 \ and J XX2 . Notice that the last three integrals are just 
derivatives of J xxx with respect to different masses. From Fig. 1 we see that there 
are two possibilities to obtain the integrals J^ d \ One possibility corresponds to a 
contraction of the lines 1 and 4 and the other one to a contraction of the lines 2 and 
3. Altogether, the representation of the integrals J^ d \q 2 ) requires us to include eight 
integrals in our minimal basis (see Appendix). 

For arbitrary q 2 the second set of recurrence relations needs a modification when 
some masses are equal to zero. If, for instance, m 2 = 0, then substituting this value 
into Eq. (|78| ) we deduce the simpler relation: 

u x (d - 2v x - 2)A 236 l + jW 2 , 3 (g 2 ) = {-Am 2 2 m 2 3 u 2 u 3 l-2 + 3 + 

+ 2m 2 2 [{q 2 - m 2 2 )(2d - 2v x - v 2 - 2is 3 - 1) - m\(2v x -v 2 - \)\v 2 2 + 
+ 2m 2 3 [(q 2 - m 2 3 )(2d - 2u x - 2v 2 - u 3 - 1) - m 2 2 {2v x -v 3 - 1)]^ 3 3 + 

-E 3 (d- vi-v-i- f 3 )W 4%„ 3 (q 2 ), (80) 

which can be used to lower the exponent of the massless propagator down to 1. There- 
fore, in this case, the number of basic integrals will be lower than in the general case. 
Each massless line reduces the number of basic type integrals by 1. 

Turning next to d- recurrence relations, we note that the reduction of d for is a 
bit simpler than that of V^ d \ The main reason is the existence of systematic recurrence 
relations which reduce simultaneously some indices and d. If two indices, say v x > 1 
and v 2 > 1, one can use the relation: 

g 2 (S 3 + 2)^ 2 l + 2 + jS 3 (? 2 ) = {q 2 {d - 2v 3 ) + m\{d -v x - 2v 2 ) 
+ m\{d — 2v x — u 2 ) — 2m 2 (d — v x — v 2 — z/ 3 ) — 2m\(q 2 — m\)u 3 Z + 



+ m\(q 2 -m\ - 3m 2 + 3m 2 - 2 + 3~)v x l 



+m 2 2 (q 2 - 3m 2 - m 2 + 3m 2 - 1" + 3> 2 2 + } J^tf) . (81) 

Relations for l+3 + J (d+2) , 2+3+J (d+2) are derivable from (|8T| ) by performing appropri- 
ate interchanges (|77|). Relations reducing the index of one line look more complicated 



-<Z 2 m 2 (£ 3 + 2)^ + l)l + l + j(S( g 2 ) + lg 2 (S 3 + 2)(£ 3 + 4)^ + J^fa 2 ) = 
{~\<l\d + 2z/ x - 2z/ 2 - 2z/ 3 ) + ig 2 [m 2 (6z/i + 6z/ 2 + 6z/ 3 - Id) 

+ m 2 (7d - 2z/ x - Au 2 - 10z/ 3 ) + m 2 (7d - 2^ - 10z/ 2 - 4z/ 3 )] 

+ m^d — 2ui — v 2 ) + m 3 (cl — 2u\ — z/ 3 ) + m\m\{d — v\ — 2i/ 2 ) 

+ m\m\{d — v\ — 2z/ 3 ) — Am\m\{d — v\ — z/ 2 — z/ 3 ) 

— m 2 [2q 4 — q 2 {2m\ + m 2 + m 2 ) + m\{m\ + m 2 ) + 3(m2 — m 2 ) 2 ]^!^ 

+ m\\q A + (g 2 - m 2 )(3m 2 - 5m 2 + 1~ - 3~) - m 4 ]z/ 2 2 + 

+ m 2 [g 4 + (q 2 - m 2 )(3m 2 - 5m 2 + 1" - 2") - m 4 ]z/ 3 3 + 

-lg 2 (S 3 + 2)(1" - 2 - 3-) - m 2 (m 2 - mjj^l+p- - 3")} 4% Ki (q 2 )- (82) 



Two additional relations follow from ( jS2f) by interchanges similar to (|77|) . 
If mi = Eq.(p2|) takes the simpler form 



g 2 (S 3 + 2)(E S + 4)^il + 4S (9 2 ) = {-? 4 (d + - 2 ^ - 2v 3 ) 
+ q 2 [m 2 2 {7d - 2^ - 4z/ 2 - 10z/ 3 ) + m 2 (7d - 2v x - 10i/ 2 - 4z/ 3 )] 
+ 2m\(d — 2v\ — v 2 ) + 2m\{d — 2v\ — z/ 3 ) — 8m2m 2 (<i — V\ — v 2 — z/ 3 ) 
+ 2m 2 [g 4 - (q 2 - m 2 )(5m 2 - 1 + 3") - m 4 ]z/ 2 2+ 
+ 2m 2 [g 4 - (q 2 - m 2 )(5m 2 - 1 + 2") - m 4 ]z/ 3 3+ 

-g 2 (S 3 + 2)(1" - 2 - 3-)} 4% U3 (q 2 ). (83) 

Again it is always possible to reduce with different shifts in d to a basic sets 
of integrals in different dimensions and then to transform all master integrals to the 



generic dimension. The relation connecting d — 2 and d dimensional integrals Jj)fl 2U3 (q 2 ) 
follow from the differential relationship given in Ref.|l|. Here we give it in terms of 
shift operators: 

JtSM 2 ) = (^1+2+ + ^ 3 l + 3+ + ^ 3 2+3+) J£l 2 Jq 2 ). (84) 

The right-hand side of this equation can be simplified by using the recurrence relations 
given previously in this subsection. As we already mentioned 4i\ and di4i\ are 
elements of our minimal set of basic integrals. For practical applications, however, it 
useful to use separate formulas for the reduction of 4i\ and di 4u- Substituting ^ = 1 
into Eq.fl8"4"D and simplifying the right-hand side we obtain the relation: 

DusAii V) = 

(d - 3) {(3d - 8)[3q 4 - 2q% - A 123 ] + nl+ + r 2 2+ + r 3 3+} j[i\(q 2 ) 

+ \mD 123 ) 4H(q 2 ) + (d 2 D l23 )J^ 2 (q 2 ) + (d 3 D 123 )4^(q% (85) 

where 

D 123 = (q 2 - ri )(q 2 - r 2 )(q 2 - r 3 )(q 2 - r 4 ) = q 8 - Ahq 6 + 2t 2 q 4 - At 3 q 2 + U, 

n = -2(q 2 - m 2 )[q 4 - 2q 2 (t 1 - Am 2 ) + A 123 ] 

r 1 = (m 1 + m 2 + m 3 ) 2 , r 2 = (m 1 + m 2 - m 3 ) 2 , 

r 3 = (mi - m 2 + m 3 ) 2 , r 4 = (m 1 — m 2 — m 3 ) 2 , 

ti = m\ + ml + m 2 = -u 123 - u 213 - u 312 , 

t 2 = 3(m\ + m\ + m 4 ) + 2(m 2 m 2 + m 2 m 2 + m 2 m 3 ) 

= 3(«i 23 W 2 i 3 + Ml 23 M 3 l 2 + M 2 1 3 M 3 1 2 ) + 2(m 2 23 + M 21 3 + u ln) 7 
i It A A A\ i 1 f 4 4 4\ , 2/ 4 4 4\ 

t 3 = m 1 (m 1 — m 2 — m 3 ) + m 2 (m 2 — m : — m 3 ) + m 3 (m 3 — m : — m 2 ) 

+ 10m 2 m 2 m 2 

= ~u 123 u 213 u 312 -u 2 123 u 213 -u 2 123 u 312 -uj 13 u 123 -uj 13 u 312 -ul 12 u U3 - u 2 312 u 213 , 



U 



A 2 



i6) 



By differentiating ( p5|) with respect to masses, and upon simplifying the right-hand 
side, we get a system of four equations connecting 4ii 2 \ ^j^iii with 4u> di4i{- 
The solution of this system is straightforward but tedious and has the form 

3q 2 (d - 3)(d - 4) (3d - 8) (3d - 10) 4ii(q 2 ) = {(d- 4)V 

- 2q 4 t 1 (d - 4)(6d - 23) + q 2 [5t 2 (15d 2 - 117d + 224) - t 2 (42d 2 - 331d + 640) 

- ~(d - 5) [* 3 (27d - 90) - t x t 2 (3d - 2) - 2t\ (5d - 26) 
+f(ml, m 2 , m 2 )l+ + /(m 2 , m 2 , m 2 )2 + + /(m 2 , m 2 , m 2 )3+} Ji ( ^ 2) (q 2 



i / 2 2 2\ T (d-2)/ 2\ , / 2 2 2\ T (d— 21/ 2 x 

+ ^(m 1 ,m 2 ,m 3 )JJ 20 '(? ) + ^(m l5 m 3 , m 2 )J^ 02 ; (g ) 



+ g(m 2 ,m 3 ,mi)4i2 2 \q 



57) 



where 



/ (m\, mj, mj) = m\(q 2 - m\) {-2{d - 4)g 4 + q 2 [u x (hd - 18) - 2Am\(2d - 7) 
-2 (id - 13) t\ + 2(9d - 31)t 2 - 24m 2 m 2 (4d - 13) - 24m 4 (2d - 7)} , (88) 



g(m\,m\,m 



1 2 
m\m\ 



A(d - 4)g 4 - 4(7d - 24)g 2 (3m 2 - 2t x \ 



3J (d - 4) 

t\ (23d - 80) + t 2 (9d - 32) - 12m 4 (d - 4) + \2m\m\(ld - 24) 



39) 



3q 2 (d - 3)(d - 4)(3d - 10)1+ J^g 2 ) = 

{-m 2 [g 4 (7d - 24) - 2g 2 ((4d - 15)*! - (d - 5)mj) 
3 5 

+ - (d - b)t\ + - (d - 3)t 2 - 2(5d - 17)m 4 - 2(13d - 45)m 2 m 2 .]l + 
+ 2m 2 (g 2 - m 2 2 )[(d - 3)(q 2 + m 2 - 5m 2 ) + (Id - 25)m 2 ]2 + 

+ 2m 2 (g 2 - m 2 3 )[(d - 3)(q 2 - 5m 2 + m 2 ) + (7d - 25)m 2 ]3+ 

- (d - 3)(d - 4)g 4 + q 2 [(7d - 30) (d - 3)ti - (Id - 3l)(3d - 10)m 2 ] 
+ \i d - 5)[(17d - 66)t 2 - (3d - 14)t 2 - 4(3d - 10)(m 4 + 5m 2 m 2 )] J jgr 2) (g 2 ) 

1 2 [(g 2 + m 2 )(7d-24) + (d-4)m 2 -(5d-18)m 2 ]m 2 m 2 42o 2) (g 2 ) 
[(g 2 + m 2 3 )(7d - 24) + (d - A)m\ - (5d - 18)m 2 ]m 2 m 2 J^ 2) (q 2 ) 
[(q 2 + m 2 + m 2 )(d - 3) + (d - 4)m 2 ]m 2 m 2 jg^V). (90) 
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Relations for 2 + j( d ) and 3 + J^ in terms of d — 2 dimensional integrals follow from (jSXf) 
by substitutions similar to ([77|). 

4.4 Two-loop bubble integrals 

Rather frequently we will encounter two-loop bubble integrals. They can be treated as 
value of 4f V2Uz (q 2 ) at q 2 = 0. It is more reasonable to consider these 
integrals separately. We found several recurrence relations most useful for practical 
applications: 

(d - 2)z/ 1 l + dS 2 , 3 (0) = {-u 123 - I" + 2- + 3-} 4 d - 2 %(0), 

(d - 2)^2+3+ J^ 3 (0) = {-2m\v x l+ + (d - 2 - 2^)} «5 3 (0), 



(d-2)(d-v 1 -v 2 -v 3 )4? V2l/3 (0) = 



- {A 123 + u l23 l~ + u 213 2~ + u 312 3~} J%-%(0), 
Ai^I+JS^O) = [ u ™( d -v x - 2v 2 ) + 2m\{v x - u 2 ) 

+ M3 i2^il + (2- - 3-) + 2m^ 2 2+(l- - 3-)} 4% U3 (0). ( 91 ) 

These relations along with those obtained from ( |9~1~D by making interchanges like ([77]) 
allow one to reduce any ^fi^(O) with integer z/j > and integer / to a combination of 

Jiii(O) and products of different one- loop tadpole integrals. Calculating integrals (]l|) 
we will encounter only two Jni(O) with different mass distributions (see Appendix). 

In the case when A123 = the recurrence relations are simpler. If all masses are 
different from zero, without loss of generality, we can assume that mf = (m 2 + m 3 ) 2 
and then by using the relation 

2m 2 m 3 {m 2 + m 3 ){d - 2v x - 2v 2 - 2v 3 - 1)1+ j£^(0) = 
{m 2 [(d -v x -v 2 - 2v 3 - l)(l+2- - 1) + [y x - v 2 + 1)1+3"] 

+m 3 [(d- v x -2v 2 - v 3 - l)(l+3-- 1) + (yi - u 3 + l)l + 2~]} jg^(O). (92) 

the integral can be reduced to a combination of products of one-loop bubble integrals. 
If mi = and m 3 = m 2 then a l so can be reduced to a product of one-loop tadpoles 
by applying the relation: 

2m\{d - 2v x -2){d- 2v x -v 2 -v 3 -\)(d- 2v x - v 2 - ^ 3 )1 + 4?U(°) = 

(d - 2v x - 2v 3 ){d - 2v x - 2v 2 ){d - Vx -v 2 - i*)jW{0). (93) 



This relation is in agreement with the exact result first obtained in Ref. P5 



4.5 One- loop integrals 

The reduction of one-loop integrals does not cause any serious problems. By applying 
the recurrence relation 



{u 612 [za1 + 2- - d + v x + 2u 2 ] + 2m\\v 2 2+\- - d + 2v x + u 2 )}G^ U2 , (94) 

along with the relation for 2 + G£p which follows from ([94] ) by interchanging v x , m x , 1+ 
and z/ 2 ,m 2 ,2 ± any G^ 2 (g 2 ) with integer z/j > can be reduced to a combination of 
Gn (<? 2 ) an d tadpoles Go^ 2 (g 2 ) and G^ (g 2 ). Thus, one-loop integrals will require 
only G11 m our se ^ °f basic integrals. The reduction of scalar integrals ([I]) will produce 



two integrals G u . One will be made from the lines 1 and 3 and the other one made 
from the lines 2 and 4 (see Appendix). To lower the space-time dimension of G^ U2 two 
relations can be used: 

2<?Vl + GX 2 = {u 126 + 1- - 2-} (95) 

2q \d - Vl -U2- 1)G { U % 2 = {A 126 + uiael" + u 21& 2~} G%- 2 \ (96) 

As in the case of the two-loop bubble integrals, it is convenient to consider one- 
loop bubble integrals separately. In principle, one can use the explicit expressions for 
one-loop tadpole integrals, however, it is more efficient to treat them on equal footing 
with the other integrals i.e. to apply recurrence relations for their evaluation. Only 
two recurrence relations are needed for all kinds of reductions: 

2mfr 1+T«> (ml) = (d - 2i*)7#>(m?) (97) 

(d - 21/OTifK) = -1m\Ti d - 2 \m\) (98) 

Five integrals T"i (mf), % = 1,..,5 have to be taken as elements of our basic set of 
integrals. Different products of one loop integrals G$ and/or which are to be 
included in the minimal set of basic integrals are given in the Appendix. 



5 Summary and conclusions 

In the present article we described a procedure for the systematic reduction of two-loop 
diagrams to a set of basic ones. We have shown how a scalar contribution of the form 
(HD, obtained from any two- loop propagator diagram, can be reduced to a sum over 30 
basic structures Ij d \q 2 ) 

30 

^ 2 ) = E^KM)4V), (99) 

j'=i 

with Rj(q 2 , {m 2 },d) being ratios of polynomials in q 2 , {m 2 } and d. All lf\q 2 ) are 
pictured in the Appendix. The choice of the elements of the basis is not unique. One 
can take some other 30 independent functions. We hasten to remark that in order to 
know the two-loop integrals of the basis one needs to deal with only three integrals 
^liiii) Win an d ^111 (^ 2 )- The other members will be obtained from those by changing 
masses or by differentiating. Two-loop bubble integrals are just the value of J\\\(q 2 ) 
at q 2 = 0. The number of basic structures strongly depends on the mass values. If 
some masses are equal to zero or there are equal masses then the number of basic 
structures substantially diminishes. For example, in the case of QED the number of 
relevant two-loop basic integrals for the photon propagator is 2, and for the fermion 
propagator it is 5. 

The evaluation of integrals from the basis is a separate problem. Two-loop bubble 
and one-loop integrals are in a sense trivial. Analytic expressions for them are widely 



known (see for example, 0], p2 |, [p6| ). As for the nontrivial two-loop integrals the 
situation is the following. Integrals with three propagators are related to the Lauricella 
functions (ITJ. As was shown in [21| integral with four propagators Vnn can be written 
in terms of multiple hypergeometric series. Similar series for -Finn with arbitrary 
masses is not yet known. In the general mass case a one-fold integral representation 
for -Finn is given in Master integrals for propagator integrals occurring in QED 
and QCD were investigated in [15], [[L7j]. Some particular cases of two-loop diagrams 
with massive particles were considered in | 23|| . We can recommended the reader for 
the general mass case to use either the one-fold integral representation |2(| or to use 
the method of Refs.i 



37 



We implemented the recurrence relations presented in this article in a FORM pi 
package. In a general mass case, required execution time for different diagrams varies 
from several seconds to several minutes. The only difficulty we encountered in calcu- 
lating diagrams with all masses different was the lengthy expressions we obtained for 
the coefficients Rj(q 2 , {mf }, d) and which are difficult to simplify by using FORM. We 
see two possibilities to solve this problem. The simplest solution will be to use other 
computer algebra system, like Maple, Reduce, etc. to simplify Rj(q 2 , {m 2 },d). The 
other solution would be to assign numerical values to the masses and then to perform 
calculations using FORM. To test the algorithm we have reproduced two-loop result 



17] for the photon self-energy. 

Several remarks concerning the possibility to extend the presented algorithm to a 
more complicated set of Feynman diagrams. Our preliminary investigation shows that 
a recurrence algorithm similar to the one described in the present article can be worked 
out also for the two-loop vertex (three-point) functions. Unfortunately, for the general 
mass case kinematical determinants which appear repeatedly in recurrence relations are 
much more cumbersome than for the propagator case. The implementation of these 
relations on the presently available computers is very problematic. For some simplified 
kinematical cases, when say, some momenta or masses are equal to zero, or several 
masses are equal, the recurrence relations are simpler and can be implemented on a 
computer. 

For integrals with several external momenta and several masses it may be reason- 
able to solve the recurrence relations numerically. It was done before in [37 using 
multiprecision package [33[] for the small momentum expansion. It was noticed that in 
this approach one needs huge arrays that frequently cause problems with a computer 
memory. However, if we do not consider small momentum expansion and work with 
unexpanded diagrams then encountered indices in integrals will be not so large and 
therefore we will not meet problems just mentioned. 
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7 Appendix 



In this appendix we present the complete set of integrals resulting from the application 
of recurrence relations for diagrams with arbitrary masses. 
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